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This observational program (in conjunction with R. B. Tully [If A, Honolulu], and J. Bland [Rice U., 
Houston]) aims to constrain the kinematic organization and dominant excitation mechanisms of ionized 
gas in dctive galaxies. More generally, we are interested in the dynamics of radiative, supersonic flows in 
the ISM. Imaging Fabry-Perot interferometers and low-noise CCDs are used for complete spatial 
coverage of the complex gas distribution in circumnuclear narrow-line regions (NLRs). Extranuclear 
emission line widths in NLRs can exceed 3000 km s' 1 , so to avoid inter-order confusion we use an etalon 
of 4000 km s' 1 free spectral range to map [N II]^A6548,6583 and Ha-C To maximize spatial resolution, 
we select nearby active systems independent of luminosity but known to possess interesting morphologies 
and/or high-velocity extranuclear ionized gas. Monochromatic images (FWHM -65 km s 1 * ) have thus far 
been obta ined in l " or better seeing at the U. Hawaii 2.2m, CFH 3.6m, and CTIO 4.0m telescopes. 

These anCst ackffij into grids of line profiles, of spectrophotometric quality, at sub-arcsecond increments 
across affifield? To handle the -20,000 - 300,000 useful spectra that arise from a typical night’s work, 
we have developed a complete analysis and reduction package for VAX and Sun image workstations. 
Reduction involves parametrization of -10 8 raw data points to a few maps (e.g. velocities of each 
kinematic subsystem, continuum-free line fluxes) containing -10 s pixels. We identify kinematic and 
structural symmetries by examining these maps and the point to point variations of the synthesized line 
profiles. The combination of monochromatic images and full spatial sampling of line profiles has allowed 
us to isolate such symmetries and has led to reliable kinematic deprojections. 


The spatially extended NLR (-6" radius) of NGC 1068 was mapped in N + with high spatial resolution 
(0."8 FWHM) at f/ 8.6 on the U. Hawaii 2.2m telescope. After subtracting the stellar continuum at each 
point, I summed the monochromatic images to generate a view equivalent to that obtained with traditional 
narrow-band interference filters. The NLR aligns with the jet axis, but correspondence to the radio 
structure is vague. The stacked images synthesize emission line profiles at sub-arcsecond spatial and -140 
km s' 1 velocity resolution. Regularities in our spectra allow decomposition of the profiles into discrete 
kinematic subsystems. These were noted by other observers in slit spectra, but flux distributions and 
kinematics can now be tracked across the entire NLR and show obvious connections with the radio 
structure (Fig. 1). Our complete coverage allows detailed study of possible interactions of the prominent 
collimated radio “jet” and associated wind with the warm ionized component of the ambient ISM. 


I find three major dynamically distinct kinematic subsystems in this NLR: 

• 75% of the total dereddened NLR [N II] flux (-12,000 [10 5 cm -3 / n e ] Mq) is contained in high- 

velocity gas (Fig. la), with approximately equal mass distributed a): to the NE, as a pair of components 

whose centroids are displaced symmetrically about systemic velocity by ±500 km s' 1 , with FWHMs 
-750 km s' 1 ; and, b) everywhere else within 4" radius, as a broad component with FWHM -1200 km 
s' 1 . Double-peaked line profiles and rapid intensity gradients are found when the jet decollimates at 

hotspots NE and SW of the nucleus. The deprojected emission distribution shows that the radio jet is 
enveloped by a thick (FWHM ~3."3 = 230 pc), center-darkened, cylinder of line emission (Fig. 2). 

The most intense gas emission is concentrated in a thick, conical shell with Gaussian centroid and 

dispersion 1/4 of the cylinder radius ( ~ 5" = 370 pc). The cone has opening angle 85° and axial 

inclination 80° to the line of sight (so its axis is inclined -40° above the disk plane in the NE 
quadrant); the jet axis lies closer to the disk plane. Intrinsic gas velocities approach 2000 km s' 1 w.r.t. 
systemic; have form r 0 - 6 with r the cylindrical radius from the jet axis; and so are inconsistent with a 
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constant velocity conical outflow. The emitting filaments are optically thick (because [O III] and [N II] 
profiles are similar out to radii of several arcseconds), with internal extinction. The surprisingly good 
agreement between modeled and observed profiles shows that profile structure is here primarily driven 
by geometry. Because no line broadening is incorporated or required in the model, the large 
extranuclear line widths in this NLR arise exclusively from geometric projection and spatial averaging 
of an asymmetric, large-scale flow, with no evidence for intrinsic “turbulence” in the outflowing 
material. If the filaments are confined by a directed, outflowing thermal galactic wind, their masses are 
> 1(T 2 M 0 , showing that they are not BLR clouds blown to observable radii. Instead it appears that 
molecular cloud cores with average densities of ~10 3 cm' 3 (as derived from CO interferometry) are 
forced by the stellar bar (prominent in near-IR images) from radii -20" into the inner NLR. There they 
are exposed to the wind, and the resulting surface instabilities strip their envelopes. The envelope 
fragments are crushed to densities of ~10 5 cm' 3 during wind acceleration, until each fragment attains 
pressure equilibrium with the shocked wind. The fragments are stable to thermal evaporation in the hot 
confining medium for the masses and velocities derived from our observations. 

• Gas near systemic velocity is concentrated in the kinematic “flare” which coincides with the NE radio 
lobe (Fig. lb); is edge-brightened, particularly along the NW edge of the lobe (leading w.r.t. galactic 
rotation); accounts for 8% of the total NLR [N II] flux (~1.2xl0 5 [10 3 cm -3 / n e ] Mq ), but negligible 
kinetic energy; and is also apparent in Ha images, so is a mass enhancement In projection this 
component is contiguous with high-velocity gas in the NE radio lobe, yet the only kinematic deviation 
from the bar-forced flow field defined by gas in the large-scale disk is localized in an expanding 
“ripple” around the base of the NE cone (Fig. Id), with kinetic energy -10 51 ergs, which appears to 
be a stand-off shock between massive, inflowing molecular clouds and their wind-accelerated fragments. 

• 9% of the total NLR [N II] flux (-1000 [10 5 crrT 3 / n e ] Mq) is emitted by a component which 
dominates within 1" of the nucleus, and extends to 2" radius along -110° P.A. (Fig. lc). Its centroid 
is blue-shifted by -250 km s' 1 from systemic, and has line width -600 km s' 1 (i.e. typical of a Seyfert 
2 galaxy). In the absence of a wind/ISM interaction we argue that this component would be the only 
manifestation of the active nucleus. Our highest-resolution images show a tight correlation of the 
brightest line emission with the nuclear triple radio-source (~0."7 extent), in agreement with speckle 
line-imaging. 

In summary, the combination of monochromatic images and full spatial sampling of line profiles allows 
kinematic decomposition of this NLR into dynamically distinct subsystems, and mapping of the flux and 
velocity distributions of each component. The global kinematic structure of the high-velocity outflow 
shows well defined, large scale symmetry. Deprojection increases kinetic energies ten-fold to ~2xl0 53 
ergs, and outflow velocities to 1900 km s' 1 . The appearance of this NLR changes dramatically as the 
velocity resolution improves from the 1000 km s' 1 (25 A) passbands traditionally employed in interference 
imaging to the 65 km s' 1 resolution attained with our current etalon. Sharply defined, bright features 
exhibit sharp kinematic boundaries, and thus resemble not static, photoionized nebulae, but rather the 
shocked termini of hydrodynamic flows. These patterns indicate that, down to the smallest observed 
scales , the spatial morphology, kinematics, and bulk of the (kinetic and turbulent) energy of the near- 
nuclear ionized gas are apparently driven by the thermalization of a directed nuclear outflow along well- 
defined shock fronts in the ambient ISM. Relativistic electrons trace only a small fraction of the outflow. 
By its spatial extension the NLR of NGC 1068 provides a unique opportunity to study the dynamics of a 
supersonic flow in an active nucleus. HST will soon image NLRs in great detail, but only with second 
generation instrumentation will it be possible to produce kinematic maps. 

We thank J. Ulvestad for the radio map and H. Thronson for the /T-band image. This phase of our 
research was supported by grants NSF 86-15974 (Hawaii) and NAS8-32902 (IAS). 
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